(Non-) geodesic motion in chameleon Brans Dicke model 
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Based on [lfj, we assume there is a non-minimal coupling between the scalar field and matter 
in Brans-Dicke model. We analyzes the motion of different matter such as, massless scalar field, 
photon, massless perfect fluid (dust), massive perfect fluid and point particle matter in this theory. 
We show that the motion of massless scalar field and photon can satisfy null geodesic motion only in 
high frequency limit. Also we find that the motion of the dust and massive perfect fluid is geodesic 
for L m — —P and is non-geodesic for L m — p. Finally, we study the motion of point particle and 
show that the motion of this kind of matter is non-geodesic. 

PACS numbers: 



I. INTRODUCTION 

Contemporary cosmology is encountered with the im- 
portant challenge of understanding the existence and na- 
ture of the dark energy component of the Universe. Anal- 
ysis of cosmological observations suggests that about %74 
of the Universe is dark energy (DE), %22 is dark mat- 
ter (DM) and the remaining part is ordinary matter [![. 
Although the nature and origin of DE are unknown for re- 
searchers until now, there are many proposals to explain 
the role of DE to explain the accelerating expansion of 
the Universe. It seems that the best one is cosmolog- 
ical constant, A, which has the equation of state (EoS) 
parameter u> = — 1, , and the second popular can- 
didate of DE model are the scalar field models with a 
dynamical equation of state. The most important dy- 
namical DE model which is called "quintessence" model 
consider the slow-roll down of a scalar field. Although, 
the quintessence scalar field cannot satisfy the local tests 
(solar system constraints), but it suggests an energy form 
with negative pressure to explain the accelerating expan- 
sion of the Universe [oT-[l2l| . 

Another suitable framework to investigate the behav- 
ior of DE is chameleon mechanism. In this mechanism 
the scalar field has non-minimal coupling with matter. 
Chameleon mechanism provides an alternative mecha- 
nism for circumventing the constraints from local tests 
of gravity. In this mechanism the scalar field acquires a 
mass whose magnitude depends on the local matter den- 
sity. Indeed this mechanism is a way to give an effective 
mass to a light scalar field via field self interaction and 
interaction between field and matter [l3|, [l4| ■ So because 
of this fact the correction of physical quantity in Newto- 
nian regime is small, i. e., the local tests are satisfied. 

Another model that has attracted much attention is 
Brans-Dicke (BD) theory. Although BD theory proved 



useful for solution of many cosmological problems, but 
it has a problem. Indeed, the Brans-Dicke parameter, 
ui, takes small value (~ 1) when standard BD model is 
used to derive the cosmic acceleration, and on the other 
hand, local constraints require that ui > 10 4 . So, some 
researchers such as, Clifton et al, [l5[ and Das et al, 
l6|, have studied another framework which scalar field 
has non-minimal coupling with both geometry and mat- 
ter. This model is called chameleon Brans-Dicke (CBD) 
model and it has predicated a value for lo which is in a 
good agreement with observational data (l7j . 

As mentioned earlier, there are a lot of attempts to ex- 
plain the positive accelerating expansion of the Universe. 
In order to do that, people have introduced various mod- 
els for example; CBD model. Actually CBD model is 
studied in detailed on large scale, solar system scale and 
it is shown that the obtained results are in a good agree- 
ment with observations fla-[23l. 
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However when a scalar field interacts with other com- 
ponents of matter (visible matter and invisible matter 
37]) through gravity or directly, this interaction may be 
produce a fifth force on the matter which may violate 
the weak equivalence principle (WEP) and creates a non- 
geodesic motion. This kind of interactions have attracted 
much attention |24j |. There are some particular mecha- 
nism for circumventing the fifth force effects. Some re- 
searchers believe that the scalar field coupled differe ntly 
to visible and invisible matter of the Universe [25H27j |. 
Therefore based on this opinion, for suppressing the ef- 
fects of fifth force, they assume that the scalar field cou- 
ples only to the invisible matter (26j . Another mechanism 
for circumventing the fifth force and then the violation of 
the WEP, is chameleon mechanism. As was mentioned 
earlier, the mass of chameleon scalar field is a function of 
local density and in the high density regions [38| the fifth 
force effects are confined to an undetectable small dis- 
tances. Therefore the violation of WEP is not observed. 
The WEP violation had become the hot topic and stud- 
ied in ditail [13J, l28| . But there are still another aspects 
of CBD model which should be studied and to the best 
of our knowledge, there is no any detailed study about 
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geodesic and non-geodesic motion for different kinds of 
matter in CBD model. 

This paper is organized as follows. In Sec. II, we con- 
sider the model, then we will obtain the equations of 
motion and conservation relation for density energy. In 
Sec. Ill, we will study the motion of massless scalar field 
and photon. In Sec. IV we consider the motion of per- 
fect fluid and point particle in this model. At last we 
summarize our work and give some discussion in Sec. V. 

Throughout this paper, the metric signature 
(— , +, +, +) and the convention 8ttG = 1 are used. 



II. GENERAL FRAMEWORK 



We begin with the chameleon Brans-Dicke action [If 



A 



d 4 x. 



<t>R - -V a <pV a <fi - V(<P) - 2/(0)L„ 



. (1) 

where R is the Ricci scalar curvature, g is determinant 
of metric, uj is the dimensionless CBD parameter, L m — 
L m {ip,g a b) is the Lagrangian of the matter and tp is the 
matter field. <\> is the CBD scalar field with a potential 
V ((/>)■ Note that the last term in the action indicates the 
interaction between the matter and an arbitrary function 
of scalar field, /(</>). 

The gravitational field equation can be derived by tak- 
ing account variation of ([1} with respect to g a b and is 
given by 

0G ab + [g ab D - V Q Vfc] 4> = f{4>)T ab + T* b , (2) 

where 

T ab = 2 ^S{^L m )0Sg ab , (3) 
is the the matter energy-momentum tensor and 

= w0</)[v a 0V 6 - l02g ab V Q </>V Q </>] - l02g ab V(<P), 

(4) 

is the scalar field energy-momentum tensor. Also by tak- 
ing the variation of (UJ with respect to tj> we have the 
equation of motion for scalar field 

(3 + 2u)D4> = f{4>)T - 4>f(<t>)L m + 4>V'(4>) - 2V(4>). (5) 

Here T is the trace of matter energy-momentum tensor 
and prime denotes derivative with respect to <f>. It is 
seen that to solve ([5]) we need an explicit form of matter 
Lagrangian, L m . 

The Bianchi identities, together with the identity 
(□V n — V a D)V c = R a b\/ b V c , imply the non-(covariant) 
conservation law 



V a T Qb = - [g ah L m + T ab ] V Q ln(/), 



(6) 



and, as expected, in the limit f(<j>) — constant, one recov- 
ers the conservation law V a T ab — 0. Since the energy- 
momentum tensor is not covariantly conserved, one may 



concludes that the motion of matter distribution charac- 
terized by a Lagrangian density L m is nongeodesic. This 
fact that the energy-momentum tensor is not divergence- 
free, can be interpreted as a violation of the so-called 
metric postulates (29j [. 



III. MATTER-SCALAR COUPLING AND 
GEODESICS OF MASSLESS MATTER FIELD 

As was mentioned earlier in the Introduction, the ex- 
plicit coupling between matter and scalar field which is 
described by the action (JTJ) can potentially lead to non- 
geodesic motion. In this section we consider the affect of 
this kind of interaction on massless particles geodesies. 



A. Massless scalar field 

Let us consider a massless scalar field ip which is de- 
scribed by the Lagrangian density 

L m = -102V.V>VV. (7) 

Using 

T a b = aL0d(v a v)v 6 v - g ab L, 

one can obtain the stress-energy tensor of the scalar field 

Tab = -V a VVfeV + l02 3ab V c ^V c ^, 19 (8) 

Substituting ([7]) and ^ in © we have 

Otl>= -V a ^V a ln(/), (9) 

One can see that for L m and L m + V c x (x is a scalar 
function), Eq. ©) is not changed just for dx,c/dV°"ijj — 
V aX,c — 0, but equation of motions for other compo- 
nents of the system, Eqs. ([5]) and ([S]), are changed in this 
case. This means that when massless scalar matter cou- 
ple with other components of the system, there is no any 
degeneracy of Lagrangian densities. 

Note that this kind of matter is known as " scalar pho- 
ton" . So, although our study is completely classic, but 
we assume this scalar photon has a wave like behavior, 
then the matter scalar field, tp, can be as a wave func- 
tion. Therefore we assume the wave function is a high 
frequency wave as 



(10) 



Here the phase of wave is a rapidly varying function of x 
and ipo is nearly constant. Therefore Eq.© becomes 

iD<f>(x) - V a $V Q $ = -iV a $V a ln(/). (11) 

Since this equation has tow real and imaginary parts, we 
have 

V a $V a $ = 0, (12) 
□$(x) = -V a $V a ln(/),Z13. (13) 
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In comparison with similar analysis in standard model 
of cosmology, Eq. (??) shows that the scalar particle is 
not transverse unless V a $ will be orthogonal to V° ln(/) 
or V a ln(/) = 0. Also by taking covariant derivative of 
Eq. (TT2|) and using V a Vh<I> = V{,V a <I>, one can obtain 



V„v Q = 0, 



(14) 



where v a = V a <I>. Note that v a = V a <I> is the gradient of 
the wave phase and in the geometric optic approximation, 
this quantity is the tangent of the worldline of the particle 
(massless scalar particle). Therefore in comparison to 
real photon, Eq. (fl"4"|) means that the motion of massless 
scalar particle in the high frequency limit is take place 
on the null geodesic. 



B. Maxwell field 

Let us consider the Maxwell field with Lagrangian den- 
sity and energy- momentum tensor 



Ln 



-I0I671-F 2 , 



T ab = -1047T 



mg ab F 2 



(15) 
(16) 



where F ab is the electromagnetism field tensor on the 
curve space-time 



F ab = \7 a A b - V b A a = d a A b - d b A a , 



(17) 



and A a is the vector potential. Substituting ([15]) and 
(HHJ) into © gives 



F ac F fc c - l04g ab F 2 



= -F ac F b c V b (lnf). (18) 



Now we consider the high frequency limit. For this goal 
we introduce the vector potential as 



A a {x) = C a e^ 



(19) 



where C a is a slowly varying vector amplitude (nearly 
constant) and &(x), the phase of wave function, is a 
rapidly varying function. So by neglecting the deriva- 
tives of vector amplitudes, C a , we have 



(v$) 2 c 2 - (V C $C* C 



2V C 



where 



A ab f{<j>) 



v a v$rc 2 - v c $c c 



0, (20) 



= 0,(21) 



A ab = C 2 V a $V fc $+ (V$) 2 C Q C h 

- (C c V c $) \c b V Q $ + C a V b $ , (22) 



and C 2 
and d22 



= C C C C , (V*) 
one can find 



V C $V C $. Using Eqs. ([20]) 



C 2 (V$) 2 = (C C V C $) 2 , (23) 
and by substituting (|2"3"1) into (|2U)) and (j2"Tj) we have 



AbV fc $ = 0, 
V b A ab = -AbV b ln(/). 



(24) 
(25) 



By setting f(<fi) — 1 one can arrive at the standard 
Maxwell equations in curved space. Geometric optics 
is valid whenever the wavelength is very short with re- 
spect to the radius of curvature of space-time, namely 
A <C £, here £ is the radius of curvature of space-time, 
and A is the reduced wave length of photon. One can 
write Eq. ([25]) as 



(26) 



\V b A ab \ ~ \A ab 0\\ » |A 6 ln(/)0£| 



So, in this case we have 



V b A ab 



0. 



(27) 



This shows that the corrections to standard optics which 
is coming from interaction between scalar field and mat- 
ter character, /(</>), complectly removed from the equa- 
tion of motion and then photons follow the null geodesies 
and are transverse. On the other hand for the case that 
X/C 1, one can not disappear the non-minimal cou- 
pling affect to the Maxwell equations and then the null 
geodesic equation of photon is modified. 



IV. MATTER-SCALAR FIELD COUPLING AND 
GEODESICS OF PERFECT FLUID MATTER 

In this section we consider a kind of matter, so called 
perfect fluid, which can be massive or massless. The 
stress-energy tensor of perfect fluid is represented by 



Tab = (p + P)U a U b + g ab P, 



(28) 



where p is the energy density and P is the pressure of the 
matter respectively, and the four velocity, U a satisfies the 
constraints U a U a = -1 and U a V b U a = 0. In [IIS] have 
been shown that, for perfect fluid that does not couple ex- 
plicitly to the other components of the system, there are 
different Lagrangian densities which are perfectly equiv- 
alent. In fact, they have shown that, by using Eq. ([3]), 
the two Lagrangian densities L mi = —P and L 
give the same stress-energy tensor as ^ 
these two different Lagrangian densities the equation of 
motions for all components of the system is similar. Also, 
since the perfect fluid laws are obtained via a kinetic the- 
ory by using microscopic models of the fluid particles and 
their interaction, namely, the perfect fluid is an averaged 
and not an exact description for matter, it is more com- 
mon to work directly with the energy-momentum ten- 
sor instead of Lagrangian density in the non-interacting 



m 2 — P 

and also for 
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model of perfect fluid. But in our model the Lagrangian 
density, L m , is explicitly appeared in equation of motion 
of scalar field, ([5]), and conservation relation, ((6]). So we 
encounter with a new situation that we have to study it 
accurately. 

We can work with stress-energy tensor of perfect fluid 
only for the case which there is no any interaction be- 
tween perfect fluid and other components of the system. 
This means that if there is a direct interaction between 
perfect fluid and other components of system, such as ge- 
ometry and scalar field, the above Lagrangian densities 
give rise to distinct theories with different predictions. 
To show this fact, let us consider a general case. 

We assume there is a minimal coupling between perfect 
fluid and scalar filed, i.e., f(<f>) — 1. In this case the La- 
grangian of perfect fluid is not appeared explicitly in the 
equation of motions of other components of the system 
and the energy-momentum tensor of matter is conserved, 
namely 

4>G ab + [ 9ab n - v a v h ] 4> = T ab + rf fc , (29) 

(3 + 2uj)D((>-(pV'(4>) +2V{(j)) = T, (30) 
V a T ab = 0, (31) 

where matter energy-momentum tensor, T ab , and the 
scalar field energy- momentum tensor, T^ b , are given by 
Eqs. ([3]) and Q. Also we suppose there are two different 
Lagrangian density L mi and L m2 which by using §3§ give 
an energy-momentum tensor like (128)1 , and L mi and L m2 
are related with together by 



A. Null geodesic of dust 

The equation of null geodesies for a model with out 
any coupling is derived from the conservation equation 
of a null dust fluid in [35[ • Therefore, in this case which 
there is a coupling between scalar field and matter we 
do the same way of derivation. If this interaction were 
to induce any corrections to the null geodesic equation, 
these has to show up in this derivation. Since dust is a 
perfect fluid without pressure, namely P — then the 
setters-energy is T ab — pU a Ub, so by using the modified 
conservation equation, (|7|). we have 



V„J7 Q = rjU a , 



where V n = U h V b and 



n 



-V„ln(//>)- V b U bl 



(33) 
(34) 



Eq. (|33|) is a geodesic equation which is non-affinely pa- 
rameterized. In fact this equation shows that the four 
velocity is transported along the path parallel to itself 
and this is the definition of a geodesic curve. So this 
means that the existence of coupling between the mat- 
ter and scalar field docs not change the equation of null 
geodesic for L m = —P. But for L m =pwc obtain 



V u U a = i 1 U a -V a Hf), 



(35) 



where clearly shows that parallel transport is no longer 
conserved and then the motion of dust particle is non- 
geodesic in this case. 



(32) 



Massive perfect fluid matter 



where x is a scalar function. Since the perfect fluid La- 
grangian is not appear in equation of motions, ()29|) and 
((30)) and conservation relation of energy, (|3T)) , then these 
two Lagrangian has not any effect on equation of motion 
of other components of the system and conservation re- 
lation of energy. Also the additional term, V a x/V~ ~9i 
in Eq. ([I]) give a surface integral term, then this term 
has no any effect on equation of motions of perfect fluid. 
Therefore these two different Lagrangian are equivalent 
for perfect fluid in a non-interacting model. 

On the other hand we assume, there is an interaction 
between perfect fluid and other components of system, 
i.e., for an arbitrary /(</>). It is obviously seen that the 
equation of motion of perfect fluid for two different La- 
grangian densities, (I32[) . are similar, moreover, by using 
Eq. (|2]) one can obtain a matter energy-momentum ten- 
sor as (l28l) . But the equation of motion of scalar field 
and also the conservation relation of energy for matter 
become equations ©and ([7]) respectively, which are dif- 
ferent with Eqs. (f2"9")l and (I3U1) and they are different for 
L mi and L m2 . This fact shows that, clearly, the two La- 
grangian density L mi and L m2 cannot be equivalent in an 
interacting system of perfect fluid with other components 
of the system. 



In this section we turn our attention to massive mat- 
ter fields and, for simplicity, consider a perfect fluid 
composed of non-relativistic or relativistic particles with 
stress-energy tensor (J25J . This stress-energy tensor is ob- 
tained from ([3]) with L m = —P. By defining a projection 
operator as h ab = g ab + U a U b , one can project equation 
© onto the direction normal to the four velocity as 



V b T ab 



[g ab L m + T ab ]V b (\nf). 



(36) 



Using (|28p. one can obtain the non-geodesic motion for 
the fluid element as 



where the extra force, f c is given by 



r = -i0(p + p ) 



(£ m + P)V 6 (ln/) + V fe P 



So by using L m = —P we have 

f = -l0{P + p)V c P. 



(37) 



(38) 



(39) 



This states the extra force is related to the coupling be- 
tween matter and scalar field and it is proportional to 
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the pressure gradient. This term is the usual term that 
appears in standard GR and encapsulates the force ex- 
erted on a fluid element due to the fluid pressure. This 
means that energy is indeed not conserved for this fluid 
does not affect the geodesic motion in GR. 

Also by substituting (|28| in to (|37|) one can obtain 

V a T ab = -(P + p)[V u (\nf)]U\ (40) 

this equation states the flow of energy only take places 
along the direction of U i.e., aligned with the fluid world- 
lines. This means that the spatial components of the 
force in the rest frame of the fluid is zero and only the 
time component of the force is nonzero. This kind of 
force cannot have any effect on the motion because, based 
on the normalization U a U a — — 1, the four acceleration 
a c is perpendicular to the four velocity U c . This states 
that the components of four-force perpendicular to the 
four velocity is zero. This is the fact which we discover 
here. Note that according to ([38| in the case of dust with 
P = 0, the extra force f° is zero and this is agree with 

On the other hand by inserting the Lagrangian density 
L m — p in Eqs. (|36|) and (|37f we have 

VaT ab = _ (P + p) [v b (ln/)], (41) 
f = -l (P + p)V c P-V c ln(/). (42) 

Equation (|41 [) states the flow of energy is not along the 
direction of U. This means that all components ( spatial 
and time components ) of the fifth force in the rest frame 
of the fluid are nonzero and then this kind of force can 
have any effect on the motion of matter. Also Eq. (|4"2")l 
shows that the fifth force does not proportional to the 
pressure gradient and it depends to the coupling between 
matter and scalar field. This means that in this case the 
motion is not geodesic. 

C. Massive matter 

In this Subsection we want to study the motion of or- 
dinary massive matter (not perfect fluid). For this goal, 
we begin with energy- momentum four- vector p<y(t). The 
density of Pg(i) is defined by (3(| 

T a0 (x, t) = po(t)5(x — x ), (43) 

where x is the general coordinate, Xo is the coordinate of 
center of particle ( the index " 0" indicates the center of 
particle). And the current of this four vector is defined 
by 

T ai (x,t)=p%(t)-l±±5(x-x ), (44) 
one can united these two definition into one as 

T ab (x,t)=p°(t)^^S(x-x ), (45) 



where x®(t) :— t. By rewriting the energy-momentum 
tensor of particle, (|45j) . in the co-moving coordinate and 
using Pq(t) — ttiqu -, we have 

T a \x) = m u a u b , (46) 

where r is the proper time, u a is four velocity vector and 

U a U a = -1. 

Moreover according to [27| , we introduce a matter La- 
grangian for a point particle with mass mo by 



L m = m Q 6(x - x )y -g ab x%x b Q , (47) 

where Xq = dx^/dr. The Lagrangian (14"7T) give the par- 
ticle equation of motion as 

x a + T a bc x c x b = 0, (48) 

This equation is a geodesies equation of motion for point 
particle. Since g a bi^x\ — g a bU a u b = — 1, one can rewrite 
(1471) in co-moving coordinate as 

Lm = mo- (49) 

Substituting Eqs. (gBJ and (ggj into Eq. © we get 

V u u b = -V 6 ln(/) - (V a u a )u b . (50) 

Eq. (HP), is geodesic equation of motion for a particle in 
CBD model which is modified with respect to conven- 
tional geodesic equation of GR. There are two terms in 
the right side of Eq. (|50l) which are the fifth force contri- 
bution from non-minimal coupling between matter and 
scalar field. Note that Vh := {g a b + u a M(,)V a is a partic- 
ular derivative in the 3-dimensional space perpendicular 
to u a . Therefore Vbln(/) is perpendicular to u a , so it 
doesn't any effect on the magnitude of velocity, namely 
this term does not any effects on the energy of particle. 
Also, since (V ' a u a )ub is aligned on the particle world- 
line, it does not any effects on parallel transport of four- 
velocity. 

V. CONCLUSION 

We have studied Brans-Dicke model which include a 
non-minimal coupling between scalar field and matter, 
so-called chameleon Brans-Dicke model. We have consid- 
ered the possible deviation of free fall trajectories from 
geodesies. We have studied the motion of massless scalar 
particles, photon, massless perfect fluid(dust), massive 
perfect fluid and finally ordinary massive particles. 

By assuming a wave like behavior for massless scalar 
field ( scalar photon), we have shown that the motion 
of scalar photon is take place on null geodesies only for 
high frequency limit, Moreover we have found that the 
electromagnetic particle (photon) is transverse and the 
motion of it is null geodesies only for the case which the 
reduced wave length of the photon be very small with 
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respect to radius of curvature of the space-time and foe 
the case A/£ « 1 the photon does not transverse and 
the motion of it is not null geodesic. 

Furthermore, we have discussed the (non)- geodesies 
motion of perfect fluid. We have found that, although 
there is a degeneracy of Lagrangian densities in the con- 
text of standard GR, but our analysis have shown that in 
the CBD model this degeneracy does not excite and also 
dust and massive perfect fluid have geodesies motion foe 
L m = — P and the motion of them is non- geodesies for 
L m = P- 

Finally, we introduced a Lagrangian and energy- 



momentum tensor for a point like particle, and we have 
shown that the motion of a massive particle is not 
geodesic in CBD theory. 
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